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We construct a macroscopic wave function that describes the Bose-Einstein condensate and 
weakly excited states, using the su(l, 1) structure of the mean-field hamiltonian, and compare this 
state with the experimental values of second and third order correlation functions. 
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The recent experimental achievement |ij of Bose-Einstein condensation (BEC) has stimulated a great revival of 
interest in the theoretical study of this phenomenon. One fascinating aspect of the Bose-Einstein Condensate is the 
, nature of coherence in a macroscopic quantum system, and in recent experiments some of the coherence properties 
of BEC have been discussed and explicitly addressed (§-§}• In this paper, we describe the hamiltonian and energy 
eigenstates within the su(l,l) mean-field picture of BEC and, based on this theory, we construct a generalised version 
of the BEC weakly excited states. We calculate some correlation functions within this theory, and compare with 
recent experimental results. 

The standard description of the Bose-Einstein condensate is by means of an order parameter field ^(x) which 
accounts locally for the physical state of the system [|| . The hamiltonian has the standard form 
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Representing the field \P(x) by its Fourier transformation leads to the second quantized form 
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H = 6 kKk + 2 V >* a t+k a t~k a P a i ( ( 2 ) 
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where Vk = V(|k|), is a momentum preserving interaction. The number operators n k = a^a k , the raising operators 
, and the lowering operators ai obey the Weyl-Heisenberg algebra commutators 

[ai,a£] = Sikl , [af ,n k ] = -6i k a^ , [ai,n k ] = 6 lk a k . (3) 

The Bogoliubov prescription is that at zero temperature the state with k = is macroscopically occupied and this 
observation allows one to treat and ao as c-numbers ([eto,a^"] ~ 0) since the corresponding number operator no, 
counting the bosons constituting the condensate, turns out to be macroscopically large. However this neglect of the 
operators and ao is not an appropriate approximation if we wish to describe phenomena in the condensate ground 



*Email: a.i.solomon@open.ac.uk 
^Email: y.feng@open.ac.uk 

^Permanent address: Dipartimento di Fisica, Politecnico di Torino, Corso Duca degli Abruzzi 24, 1-10129 Torino, Italy. Email: 
penna@athena.polito.it 



1 



states. So here we no longer adopt such an approximation and we retain the operator status of and ao in order to 
give a more consistent description of the state of the condensed system. 
Making explicit the terms depending on ao and 

in (Q) and neglecting those terms that contain three or four boson operators a^, ai (fc, I ^ 0) reduces TL to the form 

In 
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where e k = h 2 k 2 /2m (m is the fluid atom mass). 

The hamiltonian is linearized by using the mean-field approximation procedure which reduces bilinear operators 
such as AB to the linear form 

AB ~ A(B) + (A)B - (A)(B) 

based on the assumption (A — (A))(B — (B)) ~ 0. We note that a similar approach, starting instead from a bosonic 
Hubbard model and using a novel form of this linearisation procedure leads to essentially to the same su(l, 1) 
structure for the excited states, as well as an additional condensate term for the ground state. This yields the quadratic 
reduced hamiltonian 



H mf = a (n + -) + - (u a a+ 2 + a 2 ) + - E* 
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and the pair mode hamiltonian 7i^ k ' is 



y(k) _ ^-(n k + n_ fe + 1) + ^(wfca+a± fc + u* k a k a^ k ) 



with 



<r k = e k + (n )(Vo + Vfe), u k = V k (al). 
We can rewrite Tt m f in the su(l, 1) form 
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by means of the generators of the algebra su(l, 1) 
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4 0) = ^(no + 5),4 0) = ^-.^ 0) = f (6) 



and 



A k) = \{n k + n- fc + 1), Af = a+at k , A^ = a k a- k (7) 

that account for the momentum creation/destruction processes occurring in the fluid and involving the modes k and 
—k. These satisfy the usual commutation relations 

[A<*\Ato] = -2A«>, [AM,Ag>] = ±A<£\ q = 0,k,-k (8) 

It is known that within the sm(1,1) mean-field picture the energy eigenstates are expressed a direct product of 
SMfe(l, 1) coherent states 0. We therefore write the eigenstates as 

10 = |&) ®fc#o |&) (9) 

where 

|a)=exp[&4 fe) -a^ fe) ]|0) 



with th^fe = —u k /(Jk- The eigenvalues of Ji k are given by E k = \Jo\ — |itfc| 2 . 

The factor |£o) is normally absent (i.e., it is implicitly traced away) in the standard approach due to the scmiclassical 
status of Oq , do- Here it restores the condensate to its role as a dynamically active degree of freedom, that is 
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|£o} = expKo^V - C A 



where th£o = — uo/co- 
Writing for brevity 



we may express the state |£) as 



5(a)=exp[6AW-^ J ] (10) 



IO=®«Sfo)|0) (g = 0,+fc,-fc) (11) 

The operators S(^k) are similar to, but not identical with, the vacuum squeezing operators exp[i(£a +2 — £*& 2 )] familiar 
from Quantum Optics. 

The structure of the state |£) clearly exhibits the imprint of the mean-field dynamical algebra A* — (BkSUk(l, 1) 
which provides an approximate description of the dynamical processes occurring inside the system. The main unattrac- 
tive feature is the fact that 

(oo> = (flaolO = 

(arising from the two-boson character of A*), whereas the low temperature regime should be characterised by a 
nonvanishing order parameter (^(x)), or equivalently (ao)y/V (recall that |(ao) /VV\ 2 ~ N where N is the total 
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particle number inside the volume V) due to the strong depletion of the k mode states. In the state |£) we also clearly 
have (a k ) = (£|a fc |£) = 0. 

In addition, the values of the second order correlation function 

0/ 
. 2 



(a ao) 



and the third order correlation function 



ff (3) (0)= (« 

for the states |£) do not agree with the experimental results, which seem to indicate that g^ 2 \0) and g^ 3 \0) are not 
exactly equal to 1, but slightly larger than one |^,||. However it is easy to show that g^(0) — 1 and g^(0) = 1 in 
the state D(a)\0)(D state) if the mean density (no) is a large number, where D(a) = exp(aad — a*ao). 
These considerations motivate our attempt to generalise |£) to |£, z) 

\£,z) = Ifo^o) ®k^o iCfe.^fc) (12) 

by introducing the further definitions 

\^,z (s )=D{z Q )\i ) 

\tk,z k ) = D(z k )\Z k ) , (13) 

where D(z q ) = exp(z 9 a+ — z*a q ), q = 0, k, —k. 
We now describe the BEC states by |£, z) where 

\Z,z) = ® q \Z q ,z q ) = ® q D(z q )S(Z q )\0) (g = 0,±l,±2„.) (14) 

For obvious reasons, we refer to the state |£, z) as a DS state, the DS operator being similar to, but not identical 
with, that which produces a squeezed state in Quantum Optics, namely 

exp(za+ - z*a)exp[-(£a+ 2 - Ca 2 )] 

(for a single mode). 

The BEC state ( p"4"| ) involves a large number of parameters {£ kl z k } which, as is usual in mean-field theories, may 
in principle be determined by a self-consistent treatment. However, we would expect z k = for k ^ (since there is 
no condensation other than in the k = state); and if we are primarily interested in condensate properties we need 
determine only £o an d zq (4 real parameters). These may be calculated from the condensate conditions, as we now 
show. We have the following expectations: 

(£,z\a \£,z) = (£ o \D + (z )a D(z )\£ o ) = z 
(!;,z\no\Z,z) = (Zo\D + (z )n D(z )\S; o ) = \z \ 2 + sh 2 |£ | 

(£,,z\ak\£,, z ) = (6cKI6c) = z k (15) 
(£, z\a k a- k \i,z) = (£fc|a fe a_ fe |£ fe ) = z k z_ k 
(t,z\n h \£,z) = (&|n*|&) = |z fe | 2 +sh 2 |a| 
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The state |£, z) incorporates both the sm(1, 1) structure inherited from the spectrum-generating algebra approach 
to the mean field hamiltonian, as well as the nonvanishing expectation values for the operators implicit in a 
conventional (Heisenberg-Weyl) coherent state. As we shall show, a choice of the parameters for the state |£, z) state 
allows one to fit the experimental values of g^(0) and g^(0) . 

Ketterle and Miesner || pointed out that data on the condensate expansion energy, combined with spectroscopic 
scattering length measurements, can be used to give the second order correlation function g( 2 '(0) in alkali condensates. 
An experiment on a BEC of sodium by Ketterle et al yielded g^ 2 \0) = 1.25 ± 0.58, and the experiment on a 
rubidium condensate Jl(j] yielded g^ 2 \0) = 1.0 ± 0.2. In another important experiment Q, Burt et al. recently 
compared the trap loss due to three-body recombination of a rubidium condensate to that of a thermal cloud, and 
obtained 7.4 ± 2.6 for the ratio of the third order correlation function g^ 3 '(0) values in the thermal and condensed 
states. 

Although the experimental results are not inconsistent with a pure D state, at least in the case of rubidium, 
indications for sodium are that g^(0) and g^ 3 \0) are larger than 1. From the structure of the states |£, z), we can 
see that the BEC ground state is 



\£a,z o )=D(z o )S(£ o )\0) 



where 



D(zo) = exp(z a^ - Zga ) 



(16) 



),£o=r exp(i(j)) 
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with — ^5-, A_ — ^2-. We now show how to choose the parameters of the BEC ground state |£o, z a) to fit the 
experimental values of g^ 2 \0) and g^ 3 \0) cited above HQ. 

The unitary transformation of the operators ao and a,Q by D(zq) and S(£o) is given by 



S + (£ )D + (z a )aoD(z a )S(£,o) = /ia + + z 
S + (£o)-D + (zo)ao D(z )S{£o) = (J,a£ + v*a + z„ 



(18) 



where we have put 



/1 = chr, v = cxp(i0)shr 



We obtain the following mean values in the DS state: 



(n ) = M 2 + M 2 



{a+ 2 a 2 ) = 3|^| 4 + 4|z | VI 2 + M 4 + W\ 2 + »(z 2 v* + zfv) 
(a+ 3 a 3 ) = 15|^| 6 + 27\v\ 4 \z \ 2 + 9M 2 M 4 + M 6 + 9|H 4 
+ 9\v\ 2 \z Q \ 2 + 3(|z | 2 + 3M 2 M*gi/* + 4 2 ")} 



(19) 
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If we write 

h(zqV* + z*^v) = |-2 | 2 Mk| cos (^ ~ 4>z), zo = \z \exp(i<l) z ) 
Then the value g (2) (0) for the DS state is 



r 2|H 2 | \^ + \z \^W\ cos( 
(no) (n ) 



g m (0) = 1 + & + n h-imii ^npr-n^ (2Q) 



and the value g^ 3 '(0) for the -DS state is 

5 (3) (0) = l + [ 



&\v\ 2 3|^| 2 (4|^| 2 +3) 



4|^| 6 - 3(|z | 2 + 3M 2 )NVM cos(0 - 0,). 



(21) 



(™o) 3 

These results [Eg. (|20|) and Eq.(pl|)] are plotted in Figure 1 and Figure 2. From the figures we see that the 
experimental results are consistent with values of r between and 4 (r = is a pure condensate D state). 

In this note we have constructed a state for a Bose-Einstein condensate based on the su(l, 1) spectrum- generating 
algebra structure of the mean-field hamiltonian, and the Heisenberg-Weyl coherent state structure which gives non- 
vanishing boson operator expectations. It is a common feature of mean- field approximations that these give rise to 
the loss of conserved quantities (in our case loss of number conservation) and the consequent appearance of associated 
order parameters, which here are < clq > and < a,ka-k >• This is a general property of the algebraic approach 
jyj. The question of loss of number conservation is considered in some detail by Girardeau jL2|. However the 
linearisation procedure herein adopted retains the momentum-conservation properties of the original hamiltonian, as 
in superfluidity and superconductivity. The resulting DS state is similar to a squeezed coherent state, familiar in 
Quantum Optics, and will undoubtedly give rise to interesting squeezing phenomena which will be explored later. 
After completion of this work, it has been brought to our attention that a similar DS description of the condensate 
was also obtained by Navez jl3| from slightly different premises. In this note we showed that the DS state provides 
better fits to the experimental results on the correlation functions associated with the BEC state. Although it might 
be argued that the additional freedom inherent in the extra parameters associated with the DS state must give better 
fits to the experiments, it should be noted that the changes to the coherent state values (of unity) are in one direction 
only (positive) and are therefore only consistent with experimental values greater than one. 
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FIG. 1. The Second-order Correlation Function g {2> (0) for the DS state (Jz | = 50) 
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FIG. 2. The Third-order Correlation Function 5 (3) (0) for the DS state (\z Q \ = 50) 
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